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Abstract 

We investigate the properties of the path integral over reparametrizations (or the boundary value 
of the Liouville field in string theory). Discretizing the path integral, we apply the Metropolis- 
Hastings algorithm to numerical simulations of a proper (subordinator) stochastic process and find 
that typical trajectories are not Brownian but rather have discontinuities of the type of Levy's 
flights. We study a fractal structure of these trajectories and show that their Hausdorff dimen- 
sion is zero. We confirm thereby previous results on QCD scattering amplitudes by analytical 
and numerical calculations. We also perform Monte-Carlo simulations of the path integral over 
reparametrization in the effective string ansatz for a circular Wilson loop and discuss their sub- 
tleties associated with the discretization of Douglas' functional. 
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I. INTRODUCTION 



Path integrals in Quantum Field Theory are usually calculated with the Gaussian Wiener 
measure, which leads to continuous trajectories because of the presence of f x 2 in the action. 
Typical trajectories are then of the Brownian type, like those for the standard random walks. 

Another kind of the measure is associated with the path integral over reparametrizations 
or diffeomorphisms (i.e. functions with a nonnegative derivative) of a curve, which emerge 
in many interesting physical and mathematical problems, e.g. the solution [1] of the problem 
of Plateau. An important example is an integration over the boundary value of the Liouville 
field in Polyakov's formulation of string theory [2-5]. Usually, this path integral decouples [6, 
7] in the critical number of dimensions for tachyonic scalars or massless vectors, but is present 
generically for noncritical string. 

More recently the path integral over reparametrizations has appeared in an effective string 
ansatz [8] (see also Ref. [9]) for the Wilson loops W(C) in large-A QCD and was crucial in 
deriving QCD scattering amplitudes [10, 11] from these Wilson loops. The measure for the 
integration over reparametrizations was defined, roughly speaking, as 

f <T ' M 'zw--= /filer-- *'M>o, (i) 

Js(to)=s J T=TQ b \ T ) 

where r is a parameter, parametrizing the loop C, e.g. the proper length. Splitting [r , 77] 
into N equidistant intervals, a discretization of the measure (1) can be written as 

r v ms ...= *„ n r r^s f — t • • • • (2) 

Jso N ^°° 7=2^0 {Si+l-Si)J S0 (s 2 - Sl)(Sl - S ) 

where the integral goes over (N — 1) subordinated values sq < ■ ■ ■ < Sj_i < Sj < • • • < sjy = 
Sf. Discretizing the continuum formula s' = exp[— tp] and relating reparametrizations to the 
boundary value of the Liouville field v 5 ( r ) 1 by — Sj_i = exp[— ipj], we can rewrite Eq. (2) 

as 

V m s-.. = jirn^n / d^^( S/ - So -^e-^j---, (3) 

where the only restriction on the values of </?j's is given by the one-dimensional delta-function 
involved. 

1 Our definition of ip(r) is as in Rcf. [3] and differs from the standard one in string theory by a factor of 2. 
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It is clear from the representation (3) (or (2)) that there is no factor in the measure 
which makes typical trajectories to be continuous, as distinct from the Wiener measure. 
Therefore, a question arises as to what type of trajectories is essential in the path integral 
over reparametrizations, the answer to which will depend, in general, on the form of the 
integrand. 

In the present paper we investigate the properties of the path integral over reparametriza- 
tions by applying analytical and numerical methods to its discretized version (2). We refor- 
mulate the problem in the language of a proper (subordinator) stochastic process and find 
that typical trajectories are in general not Brownian but rather may have discontinuities 
of the type of Levy's flights. We study a fractal structure of these trajectories and show 
that their Hausdorff dimension varies between zero and one, depending on the form of the 
integrand. In particular, our analytical and numerical results confirm the discretization 
and heuristic consideration of the path integral over reparametrizations in QCD scattering 
amplitudes of Ref. [11] 

In Sect. II we formulate the path integral over reparametrizations as a subordinated 
stochastic process and introduce two types of regularizations of the measure: by the gamma- 
subordinator, which is convenient for exact analytical studies, and the logarithmic subor- 
dinator, which is more convenient for numerical simulations. We demonstrate that sample 
trajectories for the case when the integrand in Eq. (2) equals one, as is the case for the 
off-shell QCD scattering amplitudes [11], are expected to have zero Hausdorff dimension. 

In Sect. Ill we apply the Metropolis-Hastings algorithm to numerical simulations of a 
proper (subordinator) stochastic process and find that typical trajectories have discontinu- 
ities known as Levy's flights. Exact analytical results for the gamma-subordinator as well as 
numerical results for the logarithmic subordinator, defined in Sect. II, indicate that sample 
trajectories have a fractal structure with the Hausdorff dimension equal to zero. 

In Sect. IV we numerically investigate the asymptotic effective string ansatz [8] for the 
Wilson loops in large-iV QCD, which involves the path integral over reparametrizations, 
for a circle. Because of the presence of the Douglas integral [1] in the action, the model 
reminds ones [12-17] of one-dimensional statistical mechanics, while the difference resides in 
the measure. The expectation that the Hausdorff dimension of sample trajectories decreases 
from one at large radii of the circle to zero at small radii is supported by our numerical 
results. There are, however, subtleties of implementing numerical calculations with a step- 
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wise discretization of the Douglas integral that does not properly suppress trajectories with 
large discontinuities. We resolve this problem by using a polygon-wise discretization. 



II. INTEGRATION OVER REPARAMETRIZATIONS 

The measure in the path integral over reparametrizations (or diffeomorphisms) — the 
functions s(t) with non-negative derivatives s'(t) > — is defined as [11] 

V diS s = lim T>j$s 

N^oo 

N-l 

V^s = P(s N - sjv-i) Y[ P (si - Si-i) , s N > ■ ■ ■ > Si> Si-i > ■■■ > s , (4) 

1=1 

where s and s^ = s/ are fixed initial and final points: 

s(s ) = s , s(s f ) = Sf. (5) 

Sequences Sj with s i+1 > Sj are called subordinators in the theory of stochastic processes, 
and the sequences with fixed boundary conditions are known as subordinators that hit a 
boundary. [IS] 

The distribution P (Asj) was chosen in [11] as 

r(d i ) [Asi) 

with all equal 5t = 5 > 0. It reminds the probability density of the gamma-subordinator 2 . 
The measure for integrating over reparametrizations is approached when iV — > oo with 
N5 — > 0, while N has the meaning of the number of steps in the stochastic process. 
We shall also consider another discretization of the measure with 

P{ASi) = Hl/e t ) (A Sl + ^)' (7) 

where £j's are all the same (£j = e). In an analogy with the gamma-subordinator, we shall 
call the process caused by Eq. (7) as the logarithmic subordinator. To describe the desired 
continuum limit, we have to take the limit iV — > oo with Ne — > as will be explained below. 



2 This probability density [18] (displayed in Eq. (36) below) involves the additional exponential e ~ bAsi ; 
which cancels at the intermediate points. But it is important to provide the normalizability and an IR 
cutoff for the associated random walk problem. This factor is not essential for b(sf — sq) <C 1. 
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The consecutive AT— 1-fold integral with the measure (6) can be easily evaluated by using 
the formula 

ds, _T(5 t )T(5 t+1 ) 1 



Si-l 



We find 



S N =Sf 1 -1 

pl^-J 1 (Q) 

so 1 l 7V() J (SiV - so) 

which is an exact formula for arbitrary iV and 5 > 0. In the functional limit of iV — > oo 
with N5 we get from Eq. (9) 



[Sf - s ) 



^7 -TTM (10) 



recovering for (s/ — Sq) ~ 1 Eq. (D8) of [11]. 

The structure of Eq. (8) prompts to interpret 5 as a time step of the associated stochastic 
process, while r = NS is the corresponding time variable. Then 



dsfTv^s (11) 

J so 



describes the probability to "propagate" from s to the interval [sf,Sf + dsf] during the 
time r = N5. For small r and (s/ — s ) we can introduce the scaling variable 

z = r\n 1 (12) 

(«/ ~ So) 

which is analogous to the usual scaling variable (s/ — so) 2 / r for the Gaussian random walks, 
and rewrite the probability (11) as 

, TdS ' =dze~\ (13) 
(«/ ~ s o) 

which manifestly demonstrates scaling with 

(s f -s ) ~ e" 1 ^. (14) 

This supersedes the well- know relation (sf — so) 2 ~ r for the Brownian motion (whose 
Hausdorff dimension equals two). We conclude, therefore, that for small r and (s/ — s ) the 
Hausdorff dimension, which characterizes fractal properties of the given stochastic process, 
is zero. In Sect. Ill B we shall extend this result to (sf — sq) ~ 1. 
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We can also define the averages with respect to the measure (4), a most important of 
which is of the type 



[s K +i - s k )(sk - Sk-1, 



r s f ( (sK+i-s K )(s K -s K ^ 1 ) \ a 

Js diff ^ (sk+i-sk-i) J 



(15) 



It appears in the calculation of QCD scattering amplitudes using the principal value pre- 
scription of [11] when p(t(s)) has a discontinuity Ap K at s = s K (1 <C K <C N) but is 
constant for s < s < sk and sk < s < = Sf, so that a = a'Ap% > for Euclidean 
momenta. Using repeatedly Eq. (8), we obtain 

/f (s K+1 -s K )(s K -s K _ 1 ) Y\_ T(N5) T 2 (a + 5)T(a + 25) 

\V J /~T(a + N5) T(2a + 25) [SN S °> " Ubj 

This formula is again exact and compliments the consideration in Appendix D of [11] because 

ftfi\ NS ^ 5 r2 ( Q ) f„ „ \a 

(16) ivfM ( iV_ o) 

^ (*/ - V («/ - so) ; (s K - so) 

as is given by Eq. (D13) of [11]. 

Formulas, which are pretty much similar to Eqs. (10) and (17), can be derived also for the 
measure, given by Eqs. (4), (7), in the logarithmic approximation starting from the integral 



s i+ i ds . i n (fi+izf^i) + i n 

0* = __ 2ttl 5, (lg) 

_! - Sj + ej+i) (Si - Sj_i + Ei) [S i+1 - Si_i + £i + 

Subsequent integrations can be performed with the logarithmic accuracy using the formula 
ds- 



so 



(Sf-S + £ 2 ) (s - so + £i) 



j_ j n Ar 2 0/-s ) Jn^i- 1 ( s /~ s o) i J_ in^ 2-1 ( s -f~ s °) ln^ 1 i£iz££l 

^2 £2 £1 Ni £2 gi /iq\ 

(s f -s + £l+£ 2 ) 

which holds generically with the logarithmic accuracy for e±, £ 2 - > 0. Applying it repeatedly, 
we finally get 

fS/ #)„ 5«1 N j (20) 



^diff • 

SO 



and 



In \ (s f - s ) 

\ { (s K+1 - ) / ^ A7hTIf7>) ( 7 " o) • (21) 
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In the derivation of these formulas it was saliently assumed that N <C In i in order to 
substitute 



. V C N 1 „ / jv_! 1 



In" r. = i n ^v 1 + ^ J^v-i 1 j ; (22) 

but they may have a wider domain of applicability. Anyway, comparing Eqs. (10) and (17) 
with Eqs. (20) and (21), we conclude that 1/hi^ for the logarithmic subordinator plays, 
roughly speaking, a role similar to S for the gamma subordinator. 

With this accuracy it can be shown that the n-step transition probability 



n 



ln n_1 CnAs 



satisfies the Chapman-Kolmogorov chain condition 

P(s i+n - Sj_i,n + 1) = / dsiP(s i+n - Si,n)P(si - s;_i) (24) 

•J Si-l 

and is associated with a stable subordinator process. This can be demonstrated calculating 
the characteristic function given by the Laplace transform 

\n(eq)\ n 



J dxP(x, n) e ~ qx = ^- 



lne 



(25) 



with the logarithmic accuracy which requires qe 1. 



III. NUMERICAL SIMULATIONS 



A. The set up 

In order to simulate the path integral over reparametrizations, we apply the Metropolis- 
Hastings algorithm with the weight function (4). Note first that the distributions (6) and 
(7) cannot be normalized on the interval As G [0, +oo), which makes the usual simulations 
of free random walks impossible. That is, one cannot simulate the subordinated process 
Si by generating random As« with the probabilities (6) or (7) and adding them to obtain 
Sj. Therefore, it is essential that Sj satisfies the boundary conditions (5), which provides 
an "infrared cutoff" for the distributions (6) and (7) at large As and makes the problem 
well-posed. 

In the limits 5 — > or e — > 0, the distributions also become singular at small As. The 
measure (4) thus favors trajectories, for which all Sj are very close to each other or equal. 
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FIG. 1: Initial trajectory Si to apply the Metropolis-Hastings algorithm. 



However, since the trajectories are subject to the boundary conditions (5), some As; should 
be large. As we shall see, typical trajectories for the measure (4) consist of pieces with 
all Si almost equal and a finite number of discontinuous jumps with Asj ~ S/ — sq. Since 
the typical size of such jumps is determined only by the "infrared cutoff" (s/ — s ), they 
are somewhat similar to Levy's flights (see, e.g. Ref. [18]), which are typical for random 
walks with infinite dispersion of the step probability distribution. This is in contrast to the 
Brownian trajectories which occur for nonsingular distributions - if one simulates such an 
iV-step subordinated process with the boundary conditions (5), all steps Asj will be of order 
of (sf — sq)/N and there will be no discontinuous jumps with Asj ~ Sf — sq. It is important 
to stress that our problem of a subordinator that hits a boundary is quite different from the 
usual free random walks, where each step is not restricted to be positive. Nevertheless, for 
nonsingular P(Asj) in (4), typical trajectories will be continuous, with each step Asj ~ 1/-/V, 
in spite of the fact that our measure does not involve terms like exp (—J i; 2 ). 

In order to apply the Metropolis-Hastings algorithm, we start from the initial trajectory 
of the type depicted in Fig. 1, which is formed by small regular steps = 1/N, and generate 
Si randomly distributed between Sj_i and s i+ i. This new value is accepted if the inequality 



is satisfied for random < r < 1. The averages like in Eq. (15) are then given by arithmetic 
means over an equilibrium ensemble of trajectories. 

It is also interesting to note that since the Metropolis-Hastings algorithm only involves 




> r 



(26) 



P(s a+1 -s i (old))P(s i (old)-s i _i) 
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the ratios of the measures (4), it can be formally applied even to singular, non-normalizable 
distributions like (6) and (7) in the limits 5 — > or e — > 0. This fact should be contrasted, 
e.g., with the heatbath algorithm, which works only for normalizable distributions. However, 
in the case of a singular measure, the Metropolis-Hastings algorithm can spend arbitrarily 
large time near configurations with small As. Thus the singularities manifest themselves in 
an infinite slowdown of the algorithm. In practice, we have found that the autocorrelation 
time for the algorithm with the measure (7) grows as e~ x at small e <C (s/ — s ). 



B. Hausdorff dimension 

For 5 ~ 1 typical trajectories are similar to usual continuous trajectories of Gaussian 
random walks. Expanding the exact result (16) for large N, we obtain 

(S K+1 -S K )(S K -S K ^Y\ iV>>! N - a _ (2?) 
[Sk+1 - Sk-1) ) I 

An analogous result holds for a yet simpler average 

//, , n _ £(jVj) T{a + S) iv»i , , 

{{SK ~ SK - l] } ~ T(a + N5) T(5) ^ N " (28) 

It is easy to understand this dependence to be typical for usual continuous trajectories, when 
si — Sj_i ~ 1/N, because 

(( ( %:^r,r- J )') iy ^((^)L- 

Here r = K/N is the value of the (normalized) proper time at the point, where the derivative 
s'{t) is inserted. This scaling behavior of the gamma subordinator corresponds [19] to the 
Hausdorff dimension equal to one. 

On the contrary the result (17) in the opposite limit of TV — >■ oo, when N5 — > 0, is N 1 " 0- 
times smaller than the usual one. This holds also for the yet simpler average (28): 

;/ Na\ M5<1 (Sf — So)" / Qn x 

({sk - sk-i) ) = , (30) 

where the coefficient does not depend on a. This latter result can be understood as follows. 
In the numerical computations the average (28) is represented by 

((s K -s K „ 1 T) = /±jr^- s ^ a ) ■ ( 31 ) 

\ * = 1 / config 
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Let us assume that a given configuration has only one flight with the magnitude (s/ — so) 
that happens at an arbitrary point i. Then Eq. (30) is precisely reproduced. 

Such a configuration does not contribute, however, to the average (17) because the sum- 
mand in 



[s K +i- s K ){s K ~ s K ~i)\ a \ _ / 1 ( (sj+i - Si)(sj - gi-i) V \ 

- 1 h v - *-o J / fi { } 

I conng 



(S K +1 - S K -l) ) I \ N - 1 ~f \ - Si-i) 



would vanish. A nonvanishing contribution comes from those configurations for which a 
flight happens in two consecutive steps with the magnitudes x and 1 — x (0 < x < 1). 
An integration over x reproduces the gamma-functions in Eq. (17), while the factor of 5 is 
apparently due to the fact that the probability for a flight to happen in two steps is a factor 
of S smaller than to happen in one step. This can be demonstrated by calculating 

{(sk-sk-!)) = ^(s f -so), (33a) 

({s K +i- s K ){s K - s K -i)) = N ( 1 + N6 } ( s f - s o) • ( 33b ) 

The second formula shows the difference between the scaling regimes for N — > oo, 5 ~ 1 
and iV — > oo, N5 — > 0. These two scaling regimes correspond to the Hausdorff dimensions 
one and zero, respectively. 

A standard way [19] to calculate the Hausdorff dimension for a stochastic process is based 
on the Levy-Khintchine representation of the characteristic function 



-Tg(q) 



[ dxP(x, r) e ~ qx = e - T9[q) (34) 
Jo 



and is given by the formula 



q^oo In q 

For the standard gamma-subordinator we have 

b 



d H = hm . (35) 



P(x,r) = —-x^ 1 e- bx (36) 



and 



e -rg(q) 



= (! + §)". (37) 
so that 

lnln(l + f) t 
d H = lim ^ ^. 38 

g->oo Ing 
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This gives dn — for b ~ 1. 

For our process with discrete time r = 5N, we supersede Eq. (34) by 

^ e -q(sK+n-s K )^ — e -ng(q)/N ^gg^ 

and adopt the definition (35) of the Hausdorff dimensions, whose discrete analog explicitly 
reads as 

d H = lim — L_2 ^ ^. (40) 

q^oo In q 

It can be shown following Ref. [19] that in the limit N — > oo this definition is equivalent to 
the geometric one, obtained by covering the one- dimensional set of points Sj's by intervals. 
Generalizing Eq. (28), we obtain for (39) the confluent hypergeometric function 

e F (5n) ^ fc! r(ife + <W) i*H dn > d7V > V 

that determines at large N 

-8n jV5»l 



, -ng(q)/N 



(42) 



(l + ^)-° n "^T e-"^ for 5-1 
1 + f (e-'-l) for iV5<l 

After the substitution into the definition (40), Eq. (42) reproduces dn = 1 and dn = 0, 
correspondingly. This value of the Hausdorff dimension versus ln(l/<5) is plotted in Fig. 2 
at TV = 1000. It is clear from the figure that the desired N5 — > continuum limit is reached 
for 5 < 10~ 5 . If the value of N is smaller, this limit sets in for yet larger values of 5. 
Analogously, Eq. (25) for the characteristic function of the logarithmic subordinator in the 
logarithmic approximation results in dn — 0. 

A convenient quantity, describing the trajectories, is the spectral density 

p(x,n) = (5^(s l+n -s t -x)) (43) 

which describes how many Asj's takes the value of x. Using Eq. (42), we find 

/m . ( 5Ne- x5N (x5N) Sn - 1 /T(5n) for N5 ->■ oo 

<5 (1) (s i+n - Si - x)) = I . (44) 

[ (1 - f ) (x) + f SV> (1 - x) for 

These are two limiting cases of the exact formula 

0( x n ) = T ( N $) -ni-lQ _ \(N-n)8-l (45) 
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FIG. 2: Hausdorff dimension (40), (41) versus ln(l/<5) (left) extracted from the behavior of 
g{q) for different 5 (right) at N = 1000. The lines in the right figure corresponds to 5 = 
1, 10 _1 , 10 -2 , 10 -3 , 10~ 4 , 10 -5 from the top to the bottom. Their slope equals the Hausdorff di- 
mension. 

If additionally n8 — > oo, the density (45) has a very sharp peak at x — n/N: 



which is normalized to 1 as a consequence of the Chapman-Kolmogorov equation (24). 
C. Monte-Carlo results 

Some typical trajectories for the PDF (6) are depicted in Fig. 3 for 5 = 0.5 and S = 0.09. 
The trajectory for 5 = 0.5 looks rather normal, while the Levy flights appear for those 
at 5 = 0.09. The numerically calculated average approaches its exact value (16) after 
averaging over trajectories. For smaller S the numerical algorithm does not apparently work 
well because PDF (6) is singular at Asi = 0. 

Due to this singularity, it is very difficult to go to the required limit N5 C 1, so we switch 
to the logarithmic subordinator for which numerics works much better so we can go to the 
continuum limit for Ne 1, but analytical formulas are available only for very tiny e. 




(46) 



which is associated with the trajectory of the type in Fig. 1. 

An exact analog of Eq. (45) exists also for an arbitrary distribution: 




(47) 
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Typical trajectory for 6=0.5 
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FIG. 3: Typical trajectories Sj for PDF (6) with 6 = 0.5 (left) and 5 = 0.09 (right). 
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Typical trajectory for e=10 A -6 




FIG. 4: Typical trajectories s; for PDF (7) for e = 0.1 (left) and e = 10~ 6 (right). 

The data are presented in Fig. 4 for N = 100 and e = 0.1, 10 -6 . The values of the 
averages (31) or (32) decrease with decreasing e as is expected, but are much larger than 
ones given by Eq. (21), so this asymptote apparently does not yet set in. 

The Hausdorff dimension of the logarithmic subordinator decreases, as expected, with 
decreasing e as is illustrated by Fig. 5. 

IV. SIMULATING EFFECTIVE STRING 



We can numerically simulate the path integral over reparametrization in the effective 
string ansatz [8] 

9(C) = J V diS se~ KA ^ (48) 
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FIG. 5: Hausdorff dimension of the logarithmic subordinator versus e (left) extracted from the 
behavior of — In (exp —qAs) for different e (right) at iV = 50. 



with the Douglas integral 

A[s] 



4tt 



(ti - t 2 



(49) 



for a particular C, say, for a circle or an ellipse. Then the exponential of — ifA[s] is to 
be included in the measure. After an equilibrium set of configurations with this weight is 
generated, the averages of the type 

JV diS se- KA ^A[s] 



(A[s}) = 

can be calculated as arithmetic means. 
From the definition (48) we have 



-KA[s 



(50) 







OK 



lntf(C) 



(51) 



For a circle of the radius R, \&(C) has an expansion for large 

KB 2 = f3 



of the form 



which yields 



^(circle) = (3~ a e ( const. + 0(1/ (3)) 
A[s 



R 2 



(52) 



(53) 



(54) 



The approach to the asymptotic value tt is from above for a > 0. 
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For the sake of discretization it should be better to use a unit circle parametrization, 
changing the variable from s (—00 < s < +00) to 

a = -2arccot s (0 < a < 2vr) . (55) 

Discretizing the variable as 

0* = 2 4 (56) 



and replacing a {(pi) = Oi by 



we obtain for a circle 



N 



Si = - cot y , (57) 



A [ M ] = Jim £ sin2K (/- j0/iV] (l + ; 2 )(l + s |) ^ 58 ) 

which is a (discrete version of the) functional of Sj's over which we integrate with the same 
measure as before. The minimum of A [{sj}] is reached at 

711 

Si = - cot — , (59) 

when A min [s*] = irR 2 . 

Above Si changes from —00 to +00. For the purposes of numerical simulations it is better 
to use the variable a e [0,1], which is related to s by Eq. (57). We rewrite Eq. (49) for a 
circle as 

4W V Jo sin 2 [(0-0')/2] 
The discretized measure and the Douglas functional then read 



sin[7r(cr i - (Ti-i)\ + £ 

,.. - , 1 v-^ sin 2 [7r(o"j — cr,-)] . . 

A = vri? 2 lim — — > — 9r L ; J \ . 62 

L1 JJ n^oo (N 2 -N)*r! sin 2 [7r(i - j)/N] V ; 



The minimum of A [a] is then reached at 



o% = , ( 63 ) 



when A [a*] = irR 2 . 
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FIG. 6: Mean value of the Douglas functional (62) for a unit circle versus /3 at e 
N = 50. The solid line represents the minimal area ^4 m in[s] = tt> 
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A. Monte-Carlo simulations of subordinated process with Douglas functional 



We have performed Monte-Carlo simulations of subordinated trajectories with the mea- 
sure (61), using the Metropolis-Hastings algorithm. The average area 3 as a function of 
at e = 10~ 6 and N = 50 is plotted in Fig. 6. Typical trajectories are plotted in Fig. 7, 
and their Hausdorff dimensions are shown in Fig. 8. One can see that for small /3 the mean 
area decreases with /3, as expected. Typical trajectories become somewhat more smooth, 
and their Hausdorff dimension increases up to du ~ 0.3. However, at /? ~ 1.2 the aver- 
age area abruptly decreases below the minimal value allowed for continuous trajectories. 
Correspondingly, the jumps in the typical trajectories appear again, and their Hausdorff 
dimension decreases. 

The reason for such an anomalous behavior is that Douglas' theorem states that the 
minimal area, spanned by a given loop, is given by the absolute minimum of the Douglas 
functional in the space of all continuous reparametrizations, which are smooth. For a dis- 
cretized functional (62) this theorem is apparently not applicable. In particular, one can 
consider a degenerate trajectory with s» = if i < k and Sj = 1 if i > k for some k. For 
such trajectories the discretized Douglas functional (62) is identically equal to zero! The 



3 That is, the average value of the discretized Douglas functional (62) over an equilibrium ensemble of 
subordinators S;. 
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FIG. 7: Typical trajectories Si for the measure (61) for e = 10 6 at various (3. 
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FIG. 8: Hausdorff dimension of the subordinator with the measure (61) versus ft (left) extracted 
from the behavior of — In (exp —qAs) at e = 10~ 6 (right). 

discretized classical trajectory with Sj = i/N is still a minimum, as one can easily check 
by solving the extremum condition, but only a local one. The global minimum A[s] = 
corresponds to a trajectory with a single jump from to 1. 
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B. Path integral over reparametrizations with (improved) Douglas functional 



As we have seen in the previous Subsection, the naive discretization of the Douglas func- 
tional conflicts with Douglas' theorem. In particular, the global minimum of thus discretized 
functional equals zero rather than the minimal area spanned by the loop. The reason for 
such an anomalous behavior can be clearly identified: Douglas' theorem is valid only for con- 
tinuous reparametrizations with nonsingular first derivative, while discontinuities of the type 
of unit jumps lead to wrong values of the area. To avoid this problem, we can approximate 
the reparametrization by a continuous polygonal line rather than by discrete steps: 

a (r) =<Ti+ (<r i+ i -<Ti)(NT-i), J^< T < • (64) 

The Douglas functional can then be represented as 

= Km £ Jdn Jir ^M^MA , m 

iV^oo (iV 2 - WJ J J sm [K (ri - T 2 )\ 

where Tj = i/N and a (ri), a (r 2 ) are given by (64). The resulting integrals over the range 
[Tj,Tj+i] can then be taken numerically by dividing this interval into some number N' of 
intermediate points and applying a standard numerical integration method. The simula- 
tions would be also very efficiently accelerated if these integrals could be taken analytically. 
Unfortunately, this is not the case for the circle. Below we present the results obtained 
with thus "improved" discretization of the Douglas functional, which do not exhibit the 
anomalous behavior discussed in the previous Subsection. 

The average value of (65) is plotted in Fig. 9 versus f3 at e — 1CT 6 and N = 50. Typical 
trajectories are plotted in Fig. 10, and their Hausdorff dimensions are shown in Fig. 11. We 
can see that the mean area decreases with (3, as expected. Typical trajectories also become 
more smooth, as is indicated by a slow growth of their Hausdorff dimension. 

V. CONCLUSIONS 

Our main result is that the path integral over reparametrizations is tractable pretty much 
like the usual path integrals with the Wiener measure, while the measure is different. The 
typical trajectories have, in general, discontinuities of the type of Levy's flights and zero 
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FIG. 9: Mean value of the Douglas functional (65) for a unit circle versus /? at e 
N = 50. The solid line represents the minimal area ^4 m in[s] = ?r. 
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FIG. 10: Typical trajectories S{ for the measure (61) for e = 10 6 at various /3. 



Hausdorff dimension, as distinct from the Brownian trajectories whose Hausdorff dimension 
is equal to two for the usual path integrals. 

We have considered two types of the path integrals over reparametrizations: when the 
integrand equals one or is given by the exponential of the Douglas integral for a circle. The 
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FIG. 11: Hausdorff dimension of the subordinator with the measure (61) versus j3 (left) extracted 

from the behavior of — In (exp —qAs) at e = 10~ 6 (right). 

former emerges in calculations [10, 11] of QCD scattering amplitudes, while the latter is 
present in the effective string ansatz [8] for the asymptotically large Wilson loops. The 
Hausdorff dimension of typical trajectories for the latter is expected to decrease from one 
to zero with decreasing radius of the circle. 

Our analysis of the path integral over reparametrizations is based on representing it 
as a (subordinator) stochastic process, whose probability density has an infinite variance 
which results in several subtleties (like the central limit theorem or the law of large numbers 
are not applicable). Nevertheless, the case of the integrand equal to one can be studied 
analytically and the case of a nontrivial integrand can be studied numerically by applying 
the Metropolis-Hastings algorithm. The results justify the treatment of the path integral 
over reparametrizations in Ref. [11]. 

When the exponential of the discretized Douglas functional is included into the mea- 
sure, it makes the reparametrizations more smooth in some range of circle radii, as should 
be expected. However, at larger radii the behavior of the average value of the "naively" 
discretized Douglas functional becomes anomalous and decreases below the minimal area 
^4min[s] = vri? 2 . The reason is that Douglas' theorem is no longer applicable in the dis- 
cretized case and typical trajectories have jumps which decrease their Hausdorff dimension. 
This no longer happens, when the discretization of reparametrizing functions is polygon- 
wise rather than step-wise. Such a situation might be similar to the discretization of the 
fermionic action in lattice gauge theory, where a "naively" discretized action necessarily 
looses many properties of the continuum theory. A slow convergence of the average area to 
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its asymptote may indicate that one should think of a completely different way to integrate 
over reparametrizations in the effective string ansatz. 

Although our numerical simulations of the effective string ansatz are only for a circle, the 
method can be straightforwardly applied for other contours. It would be also interesting to 
apply numerical methods for evaluations of 2 — > 2 meson scattering amplitudes in the Regge 
kinematical regime at finite quark mass, when a spin-dependent kernel is involved [10], in 
particular, for extracting the value of the intercept of the quark-antiquark Regge trajectory. 
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